Fock-space quantum particle approach for the two-mode boson model describing BEC 

trapped in a double-well potential 
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We develop the Fock-space many-body quantum approach for large number of bosons, when 
the bosons occupy significantly only few modes. The approach is based on an analogy with the 
dynamics of a single particle of either positive or negative mass in a quantum potential, where 
the inverse number of bosons plays the role of Planck constant. As application of the method we 
consider the Bose-Einstein condensate in a double-well trap. The ground state of positive mass 
particle corresponds to the mean-field fixed point of lower energy, while that of negative mass to 
the excited fixed point. In the case of attractive BEC above the threshold for symmetry breaking, 
the ground state is a cat state and we relate this fact to the double-well shape of the quantum 
potential for the positive mass quantum particle. The tunneling energy splitting between the local 
Fock states of the cat state is shown to be extremely small for the nonlinearity parameter just above 
the critical value and exponentially decreeising with the number of atoms. In the repulsive case, 
the phase locked (7r-phase) macroscopic quantum self-trapping of BEC is related to the double-well 
structure of the potential for the negative mass quantum particle. We also analyze the running- 
phase macroscopic quantum self-trapping state and show that it is subject to quantum collapses 
and revivals. Moreover, the quantum dynamics just before the first collapse of the running phase 
shows exponential growth of dispersion of the quantum phase distribution which may explain the 
growth of the phase fluctuations seen in the experiment on the macroscopic quantum self-trapping. 

PACS numbers: 03.75.-b, 03.75.Lm, 03.75.Nt 
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I. INTRODUCTION 



The fully analytical quantum treatment of a system 
of interacting N bosons is a hard problem, especially if 
the number of particles is large, which is usually the case 
of Bose-Einstein condensate (BEC). In the limit N ^ 
1 the mean-field approach leads to the Gross-Pitaevskii 
equation for the order parameter P, [3, The mean- 
field limit may be termed "the classical limit" though the 
Gross-Pitaevskii equation describes non-classical effects, 
for instance, the quantum coherence of BEC, and reduces 
to the single-particle Schrodinger equation in the limit of 
vanishing atomic interaction. 

The relation between the full quantum and mean-field 
dynamics of BEC is a subject of intensive studies. For 
instance, manifestation of the classical bifurcation in the 
quantum spectrum of a dimer was studied in Ref. [ij. 
The mean-field dynamics of a two-mode model is shown 
to be modulated by the quantum collapses and revivals 
. In Ref. [1, 0| the Bogoliubov-Born-Green-Kirkwood- 
Yvon hierarchy for statistical operators truncated at the 
second level was used and it was pointed out that the 
many-body quantum corrections to the mean-field theory 
can be identified as decoherence. The correspondence be- 
tween the many-body quantum dynamics and mean-field 
(we will also refer to it as a quantum classical correspon- 
dence) was also studied in the phase-space by employing 
Husimi distributions Q. The classical-quantum corre- 
spondence in the context of quantum collapse and re- 
vivals in the nonlinear two and three mode boson mod- 



els was presented in [T^] , where authors described the 
BEC interband tunneling between the high-symmetry 
points of an optical lattice. It was demonstrated that 
the quantum evolution distinguishes between the stable 
and unstable classical fixed points, following the classi- 
cal dynamics about the stable points and diverging from 
it at the unstable ones. A three-mode boson model was 
treated in Ref. where a visualization of the eigen- 
states of the quantum system referring to the underlying 
classical dynamics was studied. 

The purpose of this paper is to propose an analytical 
method to study the correspondence between the classi- 
cal and quantum dynamics. We present the relation be- 
tween the classical fixed points and the eigenenergy states 
of the full quantum Hamiltonian in the limit of large num- 
ber of bosons. Our model is based on the simple fact that 
for ^ 1 the s-mode boson model is equivalent to a 
quantum mechanical dynamics of a single particle living 
in a compact (s— l)-dimensional space. We point out that 
various quantum features of the system of large number 
of bosons, not accessible within the mean- field limit, can 
be derived by the presented method. For instance, one 
can describe quantum fluctuations, BEC fragmentation 
(see, for instance, Ref. [l^]) and the cat-states formation 
in the supercritical attractive case [isl . . 

Our method resembles the classical WKB approach to 
discrete Schrodinger equation, previously developed in 
another context (see, for instance, the review [l^| and 
the references therein). The discrete WKB quantiza- 
tion method was already applied to the two-mode boson 
model [iBl- For instance, it was shown that the semi- 



2 



classical approximation reproduces the energy eigenval- 
ues even for relatively small N. 

Our approach is developed in the Fock space, where we 
expand the exact very complicated Schrodinger equation 
for the state of A^-boson system in the small parameter 
h = 1/N, the effective Planck constant, reducing it to a 
much simpler approximate equation resembling that of a 
single quantum particle in a potential. 

In general our method is applicable to s-mode bo- 
son models. Here we concentrate on the two-mode bo- 
son model, reducing it to a one-dimensional Schrodinger 
equation for the effective quantum particle. We con- 
sider the important two-mode model describing BEG in a 
double-well potential, i.e. the boson Josephson junction 
model (see, for instance, the reviews BEG in a 

double well trap is a subject of fundamental experiments 
and has many potential applications. For instance, it was 
used for direct observation of the macroscopic quantum 
tunneling and the nonlinear self-trapping il9j , considered 
as an implementation of the atomic Mach-Zehnder in- 
terferometer [iO, [2l[ , as a sensitive weak force detector 
[2^ . [2^ and for preparation of atomic number squeezed 
states for possible atomic interferometry on a chip (23 |. 

For the two-mode model we find two different 
Schrodinger equations for description of quantum dy- 
namics about the classical fixed points, corresponding to 
positive and negative mass of the effective quantum par- 
ticle. Two reduced Schrodinger equations appear due to 
the fact that a non-zero classical phase (of the stationary 
point, in our case) leads to a singularity of the limit /i — > 
and, as the result, it appears in the reduced Schrodinger 
equation. Previously, this has been accounted for by di- 
viding the turning points in the discrete WKB approach 
into the usual and "unusual" (for details see Ref. [IH). 

The paper is organized as follows. In section |TT] we 
derive the effective particle representation of the quan- 
tum boson Josephson junction model. The two resulting 
Schrodinger equations for the effective quantum particle 
are studied in section Hill The classical stationary points 
and their stability properties are studied from the quan- 
tum mechanics point of view in section IIVI In section |V] 
we derive the ground state of the quantum model. In sec- 
tion |VT] the macroscopic quantum self-trapping (MQST) 
phenomenon is related to the negative mass effective 
quantum particle. We also study the quantum dynamics 
of the running phase MQST states. Section IVIII contains 
discussion of the main results. Some of the auxiliary cal- 
culations are relegated to appendices A, B, G and D. 



II. QUANTUM PARTICLE APPROACH FOR 
EEC IN A DOUBLE- WELL TRAP 



Gonsider, for simplicity, the case of sufficiently weak 
atomic interaction, when only the degenerate energy lev- 
els of the double-well trap are occupied by BEG (we 
consider the weakly asymmetric double-well with the 
first two energy levels being quasi-degenerate, SE = 



E2 — El E3 — E2 = AE). Under this condition one 
arrives at the reduced Schrodinger equation (see also ap- 
pendix S and Refs. [li,[2|) 



-ff = I -^{a\a2 + 4ai) + + -^{ni + n^) 



(1) 



N 



where \ip) is the system state, rij — OjUj (the bo- 
son field operator is expanded over the localized states 
^ = (/3i(x)ai 4- ipji{x)a2) and the time is measured in 
the tunneling time units: r = t/T. The parameters of 
the model are defined as 



2n _ 26V 
~6E' ^^~5E 



9N I ,3 4 



j d^x [ipLV{^)ipL - ^RV{^)ipR] . (2) 



sv 



where g = Anhas/m with as being the s-wave scatter- 
ing length, V(x.) is the double-well potential (centered at 
zero, for below), the localized states ipl{^) and V'fl(x) are 
given by the normalized sum and difference of the ground 
state and the first excited state of the "symmetrized" 
double-weU potential 14 (x) = [V{x) + V{-x)]/2 (due to 
the symmetry the integrals of the fourth power are equal; 
for more details, see Appendix in Ref. [26|). 

The dimensionless parameter SV can be of the same 
order as SE, thus e can be arbitrary. Note that 7 is finite 
(not small) and the applicability condition that the inter- 
action energy is much less than the average trap energy 
spacing, i.e. 7 <C AE/SE, is satisfied due to the large 
r.h.s.. 

Subtracting the conserved quantity 2^("'i+"-2)^ from 
the Hamiltonian in Eq. ([1]) we obtain the standard rep- 
resentation of the Hamiltonian for the boson Josephson 
junction (see, for instance, Refs. @, Ul]). Our interac- 
tion parameter is related to the standard parameters as 
follows 7 = EcN'^/Ej. 

We have divided the Schrodinger equation ^ hy N 
and defined the parameters by extracting the explicit 
iV-factors reflecting the "order" of the respective boson 
operators (e.g. a\ai ~ N). This is used in the effec- 
tive quantum particle representation of Eq. ([T]) (see Eq. 
([5])). It is important to observe for the following (and 
the validity of the explicit iV-factors) that the physical 
parameters of the model are TV-independent in the ther- 
modynamic limit, which is — > 00 at a constant density 
where f2 is the volume. Indeed, T and e are evidently 
constant, while the nonlinear parameter 7 '--^ since 

On the other hand, the factor N^^ at the time deriva- 
tive on the l.h.s. of Eq. ^ plays the role of an effective 
Planck constant and the thermodynamic limit iV — > 00 
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now acquires another meaning, namely it is the semiclas- 
sical hmit of Eq. ([T|). To make this clearer, let us rewrite 
Eq. ([1]) in the explicit form using the Fock basis 



N 



IV) =^C,.|fc,iV-fc), 



(3) 



fe=0 



where the occupation number corresponding to the left 
well is denoted by k. We get 



i d 



{bk-iCk-i + bkCk+i) +akCk, (4) 



with 



bk = [{k + 1){N - k)]-^ , 

ttk = -/N-^ [k^ + iN- kf] + sN'^k. 

Introducing a continuous variable x = k/N , the effective 
Planck constant h — 1/N and a wave-function \E'(x) = 
VlV Cfe we obtain Eq. ^ in yet another form 

ihdr^{x) = - {e-'Pbh{x) + bh{x)e'^] ^{x) + a{x)^{x), 

(5) 

with p = —ihdx and 

bh{x) = [{x + h){l^x)Y/'^, a{x) = -^[x^ + {l~xf]+ex. 

Here x G (0,1) and p are canonical variables with the 
usual commutator [p^x] = —ih [35]. Note that the inner 

product for TV » 1 reads {^i\'^2) = dx^l{x)'^2ix). 

Eq. (O is written in the ^-representation form. It can 
be put in the (/)-representation by putting p ~f (p and 
X — 1/2 ^ ihd^ (the shifted variable is more convenient). 
In a more rigorous way, one can use the transformation 
of the state vector Ck to the Fourier space, given in our 
case by the discrete Fourier transform, 



1 ^ 



-'''^'Ci, Ci 



1 ^ 



(6) 

where (pi — 2ttI/{N +1). Using periodicity of the expo- 
nent we have for ^ 1 



where {...} denote the anti-commutator. The Hamilto- 
nian ([5]) reduces to the usual non-rigid pendulum form 
(27| if one neglects the variation of bii{x) with respect to 

X. 

The limit of large N of Eq. ([5]) can be considered by 
the WKB approach for the discrete Schrodinger equa- 
tion [iBl . We write the solution of Eq. ([5]) in the form 
\I' = exp{iS{x,T,h)/h} where S(x,T,h) is understood 
as a series S = S^°\x,t) + hS^^\x,T) + OiK^) and the 
terms are assumed to be differentiable functions. In the 
lowest order we obtain from Eq. (O the Hamilton- Jacobi 
equation for the classical action S''^*'^ (a;, r): 



Sf^ =a{x)-2b^{x)cos(^^\ 



(9) 



where ba{x) = \J x{\ — x). The classical Hamiltonian 
obtained from Eq. ^ reads 



H = a{x) — 2bo{x) cos (j) 



(10) 



with the canonical variables x and (p e (0, 27r): {0, x} = 1 
(cf. with Ref. [27], where z ~ 1 — 2x). The first-order 

approximation reads S*^^^ = | ln[6o(a;) sin(5i°^)]. 

Let us, however, take a slightly different approach. The 
most important features of the classical dynamics are the 
stationary points and their stability properties. About a 
stationary point (say, x^ and pg) the classical action can 
be expanded as follows 

S^°Hx, t) = -E^-'\ + - xs) + 0[{x - xsfl (11) 



where we have used that S*'"^ (x, t 



and 



-{Xs) = Ps 



-E('=^'>t + S^'*\x) 
We are interested in the limit of 



large N. The wave-function localized about some Xg- 
point has a singular derivative in the limit h 0: ^ ~ 
ih~^(l)s'^{x,T) if (ps 7^ 0. The singularity is due to the 
classical phase in the wave-function. We then account 
for the phase explicitly, introducing the transformation 
(at some x^-point) 



^{x,t) 



JN<I,{t){x-Xs) 



ip{x, 



(12) 



^(a;) J d^e*^"^^^ (</.), = AT J dx e"'^'^^*(a;), 

-TT 

(7) 

where ^{pi) = {N+l)Ci. Thus, for TV > 1 the Hamiho- 
nian for boson Josephson junction in the ^-representation 
reads (save an inessential constant term) 

H = -2-fh^dl~{b,,{l/2 + ihd^),cos{<P)} 

+iehd^ ~ i[bh.{l /2 + ihd^), sin (j)] (8) 



before taking the limit of large N. The phase 0(t) sat- 
isfies the evolution equation 4> = dH{xs,4>) / dx (by ex- 
panding the Schrodinger equation in x — Xg and taking 
the limit /i — > 0), hence it is nothing but the classical 
phase. 

To obtain a reduced Schrodinger equation we substi- 
tute the representation (fT^ into Eq. (O and expand the 
result into series in p keeping the terms up to OijP') (since 
the derivative is regular the expansion is actually in h). 
We have 



4 



dx 



'"^fbhix) + e'Hh{x)f] + 0{f) Uix) 



(13) 



r 



One can make one more reduction in Eq. (fT3|) by using Using Eqs. and (fTil) in the Schrodinger equation ([51) 
that we obtain 



1 



-^^sin^(:^^ + p:^l+-cos(^)^g.r(i4) 



ihdril^ = V/i(x, 0)V' + cos 4>pbh{x)p + sin </> + j56h(a;)) - sin0 ("^P + ) ^ + O(p^) 



(15) 



r 



where the potential Vh reads 

Vh = a{x)~2bh{x) cos <j) + h cos (f> 



dx 



4 dx^ 



(16) 



Reducing further, one can neglect the terms of order 0{h) 
in the potential Vh, since it contains the /i-independcnt 
term, obtaining nothing but the classical Hamiltonian 
(Uni) Vo(2:) = H = a{x) - 2boix) cos (j). 

The analogy with the quantum problem of a parti- 
cle in a potential becomes evident from the Schrodinger 
equation in the form of Eq. (fT5|) . The stationary states 
ij{x, t) = e~''^'^^'^^pEix) correspond to the extrema of the 
quantum potential. The extremal values of the classical 
phase (which is a parameter) given by ^''^"g^''^-' = co- 
incide with the classical values at the stationary points: 
(j) = 01 (j) = TT . These correspond to positive or negative 
mass of the effective quantum particle. 



III. QUANTUM MODEL ABOUT THE 
MEAN-FIELD STATIONARY POINTS 



The Schrodinger equation (|T5|) in the two cases (/)+ = 
and 0_ = TT reduces to 

ihdri^ = {±P\/x{l - X)p + V± [x] } V^, (17) 

where the upper/lower sign corresponds to <f>+/4>- and 
V±{x) = 7[x^ + (1 — x)'^] + ex =F 2-y/a;(l — x). For further 
analysis it is convenient to introduce a new variable z = 
1 — 2x, the atomic population difference divided by the 



total number of atoms. We get 



dz 



z^^ + V±{z)}^b, (18) 



where the potential simplifies to 
1 



V±{z)^-hz' 



ez 



(19) 



Note that the potential V±(z) crucially depends on the 
nonlinearity parameter 7. 

Taking into account the symmetry of Eq. (llSp with 
respect to replacing ± we have the following mapping 
between the repulsive and attractive BEG cases: 



7attr 



-7ro 



Zattv 



i'eittT{z,T;-/,(j}±) = ■(/'rep(-^,'r;-7,(/)zp), (20) 

i.e., for instance, an attractive BEG with positive mass 
^ = is equivalent to a repulsive BEG with negative 
mass (j) = TT evolving backwards in time in the inverted 
space z — > — z. Note that it is sufficient to consider also 
e > 0, since £ < can be simulated by replacing z — > —z. 

The equivalence ([201) . however, is useful only in the 
general analytical analysis, the attractive BEG is funda- 
mentally different from the repulsive BEG. As we show 
below, in both cases the Schrodinger equation p8|) with 
the positive mass gives the ground state of BEG and only 
in the attractive case V+(z) has the double- well form (see 
Fig. [Il^a) below). 

The extrema of the potential solve the equation 



7± 



1 



e 

Yz' 



(21) 



5 



Consider 7 > 0, e > 0. In the case of positive mass 
((/) = 0) there is just one minimum. For the case of nega- 
tive mass (0 = tt) there is a critical below which there 
is just one solution which is a minimum ofV = —V-{z), 
whereas for 7 > 7cr there are three solutions, two corre- 
sponding to local minima and one to a local maximum, 
see Fig. [Tfa). Eq. (PTjl is easily solved graphically by 
plotting the l.h.s. and the r.h.s. as functions of z, see 
Fig. [TJb). Hence, in the case of ^ = (the upper sign) 
we get that the positive-mass effective quantum particle 
moves in a single well potential V ~ V+{z) and there is 
only one solution to Eq. (PT|) . In the case of cj) — tt the 
negative-mass particle sees either one well or two well 
potential V = — V-{z). 
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FIG. 1: (Color online) The potential for the effective quantum 
particle, panel (a), (here V = —'V^(z) for the repulsive BEC 
corresponding to negative mass, </> = tt, or 1^ = V+(— 2;) for 
the attractive BEC for </> = 0, positive mass) and the graph- 
ical solution of Eq. (I2ip . panel (b) (again for = tt for the 
repulsive BEC, or for </> = for the attractive BEC, in the 
last case the axes are inverted). Here e = 0.3 which gives 
|7cr| = f.452f. 



The critical value of the nonlinearity coefficient can be 
found by equating the derivatives of the two sides of Eq. 
([2T|) . which gives an additional equation. We get 



(22) 



while the corresponding solution, i.e. the point of inflec- 
tion, reads Zcr = [1 + (2/e)^/'']^-^/^. For 7 > two pos- 
itive solutions bifurcate from it: z\ < Zcr < ^2, where z\ 
is the local maximum of the inverted potential — VL(z). 
The local maximum solution is an analog of the z = 
solution for e ^ 0. 

One can show (see Appendix |B]) that in the case of 
(j) = TT the negative local minimum of the inverted quan- 
tum potential is also the absolute minimum (of this po- 
tential, i.e. —V-{z)) for all 7 > 7cr and e > 0, whereas 
the negative z correspond to higher atomic population of 
the trap well with the higher zero-point energy (see Eq. 
P])). This fact, however, is easily explained within our 
approach by the negative mass of the effective quantum 
particle. 



Finally, all the results obtained for the repulsive BEC 
are transferred to the attractive BEC case by using the 
equivalence given in Eq. ([^0]) . For instance the cases of 
(j)± are interchanged and the double- well potential of Fig. 
[Ija) in the case of attractive BEC appears for (f> — (i.e. 



now V = V+{—z)) and 7 < —7c 



f 



(§) 



2/3 



3/2 



IV. THE MEAN-FIELD STATIONARY POINTS 
AND THEIR STABILITY EXPLAINED 

Before we proceed with the analysis of the bound state 
wave- functions, let us explain from the quantum mechan- 
ical point of view the appearance and stability properties 
of the stationary points of the classical Hamiltonian (see 
also Ref. [131 )■ In the new variables we have 



(23) 



where one must remember that {(f), z} — —2. The 
stationary points of the classical Hamilton equations, 
z = —2^, (j) = 2^, correspond to the extrema of the 

quantum potential V{z), i.e. cos 0s — ±f, whereas Zg is 
defined from Eq. (HJ). Consider the case 7 > (repul- 
sive BEC) and e > 0. The stability properties are defined 
by the local Hamiltonian, see also Fig. [21 Introduce the 

local phase space co-ordinates cf) = (jj^ + ^p and z = Zs + C- 



y<r 



Y>Y„ 




FIG. 2: (Color online) The phase portrait of the classical 
Hamiltonian Eq. (|23p corresponding to a repulsive BEC in a 
double-well trap. The parameters are as in Fig. [T] 



I. In the case of = we have up to the second-order 
terms in the local variables 



H+=V+{zs) + VT^s^ 



7+(l-^') 



'3/2 



2 



(24) 

hence the classical stationary point (pg = Q and z = Zg, 
solving Eq. pT|) for the upper sign, is elliptic and stable. 
This stationary point, in fact, corresponds to the absolute 
minimum of the quantum potential V+(z). 

2. In the case oi (pg = tt we obtain up to the second- 



6 



order terms (by simply changing the sign at cos ( 



V * 2 



-3/2 



■7 



V-{zs) 

(25) 

First, consider the sohition to Eq. (pij) such that < 0. 
Using Eq. pT|) in the following estimate 



1 



(1 



z2)3/2 



< 



7 



1 - Z2 (1 _ ^2)3/2 



2z,(l-z2) 



<0, 



we conclude that the negative stationary point < 
for 0s = TT is always elliptic and hence stable. It actu- 
ally corresponds to the absolute minimum of —V^{z) (see 
Appendix [B)) . 

Now consider the positive stationary points z^^"* < 
Zcr < z"^\ which appear for 7 > 7c 



For z's'^ we have 



1 



(1 - 



3/2 



>7- 



(1 



y2 \3/2 



7 - 7cr > 0, 



where we have used that 7cr = (1 — z^^)""^^^- Thus this 

stationary point is hyperbolic and hence unstable (recall 

that it corresponds to the local maximum of — F_(z)). 

(2) 

The other stationary point Zs > Zcr is a local minimum 
of the quantum potential —V-{z). Let us show that it 
is always elliptic, i.e. stable. Indeed, using Eq. ([2T|) we 

(2) 

obtain the estimate (z = Zg ) 

1 / 1 



7 - 



(l-z2)3/2 
1 / £ 

Yz 



7 - 



1 



7z 



1 - z 



1 - z2 
'"2^3 



1-Z2 

< 0,(26) 



1-Z2 

> Zcr, the expression 7cr = 



(2) 

where we have used that Zg 
e/{2z^j.) and that 7 > 7cr. 

According to Eq. (jSO]) , the classical action for Eq. (flSl 



satisfies sI^I^{z,t) — — 5rcp(— z,t) where 7 



(0), 



-7 and 



thus Eqs. 



and can be transferred to 



an attractive BEC by the substitution: —Ti± 

-v± 



Tizp and 



In conclusion of this section, the local minima of the 
quantum potential (or the inverted potential, in the case 
of negative mass) correspond to the elliptic stationary 
points of the classical dynamics and the local maximum 
to a hyperbolic stationary point, which is a physically 
clear result. 



V. THE GROUND STATE OF BEC IN A 
DOUBLE- WELL TRAP 

To find the ground state of BEC in the double-well 
trap one has to analyze the energy of the bound states 
localized at the extrema of the potentials V±(z). For 



N ^ 1 one needs to compare just the zero-point energies 
at the extremal points (which is the classical energy of the 
stationary point, Eqs. ([M)) and (^5]) ). The result is that 
for both attractive and repulsive BEC in a double-well 
trap the ground state is given by that of the effective 
quantum particle with the positive mass, i.e. with the 
classical phase = (see appendix [C] for details) . 

The direct link between the stable classical station- 
ary points and the nature of the corresponding quantum 
states discussed in section IIVI allows one to get the local 
approximation for the quantum bound states by directly 
quantizing the local Hamiltonian, Eqs. and ([^5]) . by 
replacing if — > 2ihd(^, since [(p,C] = 2i/i (remembering 
that in the case of attractive BEC Eq. p4)) corresponds 
to = TT, while Eq. ^ to = 0). We shghtly cor- 
rect this scheme in the negative mass case, where the 
wave function '^{x, t) of BEC in a double- well trap has a 
non-trivial phase according to Eq. (fT2l) . 

Writing the local classical Hamiltonian in Eqs. ((24|) - 
(Uni) as 

H±[C.V) ^V±{zg) + [a±'p' + b±e)/'^. (27) 

we obtain an approximation to the bound state of BEC 
corresponding to a stable stationary point of the classical 
dynamics and its energy as follows: 



1 



1/4 1/2 



exp 



.0± {X-Xs) 
1——X 

h 



2a2 



± 



e± = v±{zs) + o{h), 



(28) 



is the Gaussian width and we have 



where a| = ^^^f^ 

used that Az — 2Aa;. The energy spacing for the few 
lower levels reads A£± — hU±, where n,± — 2y/a±b± 
is the classical frequency. In the case of negative mass 
the energy levels of a local quantum Hamiltonian are de- 
scending. 



A. The ground state of repulsive BEC 

For 7 > we obtain the single solution of Eq. (pij) 
corresponding to the positive mass case in the form of a 
series in e: 



= 2(1 + 7)- 



1 - 



8(1 + 7)' 



(29) 



which is valid for a weakly asymmetric potential V+(z), 
i.e. £ <C 1 -f- 7. The average atomic population differ- 
ence between the two wells reads {ni — rt2) = ZsA^. By 
the substitution 7 — > —7 and z^ — > —Zg Eq. (j29p also 
gives the position of the minimum of the inverted poten- 
tial — V_(z) in the subcritical case 7 < 7cr- 

For a weakly asymmetric potential V+(z) one can also 
calculate the atom number fluctuations using the local 
approximation for the ground state given by Eq. 



7 



2^2 



£'(4 + 7) 
16(1 + 7)3 



where we have used that 



2N 



(1 



\l/2 



7+(l-z2)-3/2 



1/2 



(30) 



(31) 



and the series expression ((29)) (see the details in Appendix 

El). 

The resuh given by Eqs. ^ and ^ can be 
compared to the non-interacting case, which is exactly 
solvable (here we consider the case e = 0). Indeed, 
the tunneling term of the quantum Hamiltonian ([T]) 
can be diagonalized by the canonical transformation 
c± — (oi ± a2)l\f2 with the effect H — — (a|a2 + 
a}^^-^)!^ — (clc_ — c^c_|_)/A^. The eigenstates are given 
by (c^^)™(cl)^-''-|0,0), where ai|0,0) = a2|0,0) = 0. 
The ground state corresponds to m 



N 



1$) = 



1 



N 



|o,o). 



(32) 



For large N, approximating the factorial, we get the co- 
herent state in the Fock basis as 



<f(z) 



N{ 
2N 



1 



exp 



N 

T 



+ 



6 



(33) 



Eqs. ((28)) - ((30)) describe all three known regimes of 
repulsive BEC tunneling in the double-well trap (see for 
instance, Ref. [l3,[lB|): (1) Rabi regime, when the co- 
herence is very high and the atom number fluctuations 
are large (essentially the interaction free regime), 7 <C 1; 

(2) Josephson regime, when the coherence is high and the 
atom number fluctuations are small, 1 <C 7 ^ N'^ and 

(3) Fock regime, when the coherence is low and the atom 
number fluctuations vanish, 7 3> iV^. 



B. The ground state of attractive BEC 

Subcritical 7. In the case I7I < 7cr the ground state 
of attractive BEC is essentially the same as that of the 
repulsive BEC and is given by the general result ((^5| for 
the positive mass case (i.e. the upper sign). The average 
population difference is again given by the series solution 
(|29p and the number fluctuations are given by the same 
formula ((50)) as in the repulsive case. However, since now 
7 < 0, there is an essential difference: even moderate at- 
tractive interactions strongly enhance the atom number 
fluctuations (which fact leads to absence of the Josephson 
regime for attractive BEC). 



For 7 approaching (from above) the critical value — 7cr 
the fluctuations given by the local oscillator approxima- 
tion ((^5)) diverge, which is easily seen from Eq. ((50)) 
applied to the case e = 0. In the general case the result 
is similar and follows from the exact expression for cr^ 
with the use of the exact value for the point of inflection 
Zs (given below Eq. (1^^ of section IIIip and the expres- 
sion 7cr = (1 — ^cr)"^^' The divergence results from 
the oscillator approximation of the wave- function \['(a;) 
((^5)) which breaks down at the critical value of 7, since 
the potential V+(z) for 7 = — 7cr is approximated about 
the stationary point z — Zg hy the fourth order power in 
z — Zg (we give the case e = and = for simplicity): 



V+iz) 



y2 - 



Z 

16 



(34) 



Hence, the usual oscillator approximation should be re- 
placed about the critical 7 by the ground state in the 
non-harmonic potential, not known in the exact form. 

Supercritical 7. In the general case of I7I > 7cr one has 
to solve the stationary Schrodinger equation ((T^ to ob- 
tain the ground state. However, for 7 such that the po- 
tential V+{z) is a double-well with well-separated wells 
(the form of this potential is essentially determined by 
the atomic interaction parameter), the oscillator approx- 
imation is still valid locally (i.e. about the two minima 
of the potential V+{z), see Fig. [Ij and the ground state 
can be now approximated by a combination of the lo- 
cal eigenfunctions of the form ((28| due to the quantum 
tunneling between the wells of V-f (z). 

The case £ = is the most interesting. The two local 
minima, solutions of Eq. ((21)) . read Zs*-* = ±-\/l — 7^^. 
The validity of the local approximation by oscillator 
eigenfunctions in the two wells, i.e. the condition of well- 
separated wells of the potential V+{z) for 7 < —1, is 
defined by how much the width of the local oscillator 
eigenfunctions is smaller than the distance between the 
two wells. We have zb. — zl — 2-\/1 — 7"^ and using Eq. 

{{z-ZL?) = mx-XL)')^N-^\j\-\j^~ir^/\ 



Therefore, the applicabihty condition ((z — ZLf) <C 
(zfl - zl)2/4 reads 



I7l(l7|-|7r^)^» 



(35) 



where on the l.h.s. we have a monotonically growing 
function of |7|. The condition is not too restrictive and 
already for I7I slightly above the critical value I7I = 1 it 
is satisfied, e.g. setting I7I = 1 -I- (5 with (5 ^ 1 in Eq. 
((55)) we obtain 5 > which for N = 1000 gives 

S ^ 0.01. Hence already I7I = 1 + (5 > 1.1 satisfies the 
condition ([55]) for iV = 1000. 

To obtain the ground state we need to evaluate the tun- 
neling rate, i.e. the matrix element (^^jiJl^'fl), where 
for the localized eigenfunctions in the left and right well 
of the double well V+{z) one can use the local approxi- 
mation ((^5)) and H can be taken from Eq. ((15)) . We have 
for the localized states in the left and right wells (here 
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1 



2 |„|-U,,2 l)-l/2_ 

is an exponentially small overlap between 
1 



where tr^ = 4f7^ = ;^|7| ^(7^ 



(*L|*fl) 



1 



dz exp{- 



(2 - ^l)^ 



1)3/2 



where we have used that {z — + (-^^ ^ '^nf = 

2(^2 + z|) and the normalization of the function '^{z) = 
7r^^/''cr^^/2 exp{— ^}. Using the integration by parts 
we get 



- zl){z - zr)VI-z^\^r) 



To estimate the products we use (z — zl){z — zn) — _ 
1 + 7~2 and set z = where the overlap of the wave- 
functions l,r{z) is maximal, with the effect [ST] 



1 



(38) 



It is natural to define the "tunneling coefRcient" k for the 
effective quantum particle as follows (cf. with Eq. ^) 



(72 - 1)2 
1 + ^ '- 



(39) 



exp 



TV 
"2H 



(7^ - 1) 



3/2 



The ground and first excited states in the supercritical 
double-well potential V+{z) for 7 < read 



\G) 



1 



\^r)), \e) 



1 
71 



{\^l)-\^r)), (40) 



where ^l,_r(2;) are given by Eq. ([55]) . The theoretical 
ground state is a very good approximation of the numer- 
ical diagonalization of the quantum Hamiltonian ([1]), see 
Fig.H 

The degenerate level energy splitting, given by SE = 
2k, is a qualitatively good result, see Fig. ID Note that 
the A^-dependence of the energy duplet splitting is mani- 
festly exponential (however, the local approximation used 
for the wave functions cannot capture the right coeffi- 
cient in the exponent). The energy splitting can be also 
estimated by using the perturbation theory, but only for 
small values of N (since the small parameter is equivalent 
to N/-f see, for details, Ref. [28*1). 

The atom number fluctuations in the ground state (|40|) 
read 



(Anl) 



7V2 



(G|z2|G) 



7V2 



(I-7-'). (41) 




FIG. 3: (Color online) The ground state of attractive BEC 
in the supercritical case: theoretical (given by Eqs. (|36p and 
(|40p . solid line) vs. numerical diagonalization of the Hamil- 
tonian (dashed line). The thick grey line shows schematically 
the effective potential V+{z). Here A'^ — 200, 7 = —1.2 and 
e = 0. 
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FIG. 4: (Color online) The energy splitting due to tunneling 
in the effective potential well vs. number of atoms A'^: theo- 
retical (solid line) and numerical (dotted line). Here 7 — —1.2 
and e = 0. 



For large I7I above the critical value the ground state 
of an attractive BEC in the symmetric double-well trap 
is the Schrodinger cat-state (see also Refs. [l3,[l3l)- Pre- 
cisely, for I7I 3> Vn we have 



1 



|G) = ^(|A^,0) 



|0,7V)). 



(42) 



Indeed, in the supercritical case '^l,r{z), given by 
Eq. (1551) . has the width in the Fock space given 



by 7VV((z-z,)2)/2 
%/iV/(2|7|) < 1 for I7I > 
is a Fock state. 



I 



^ ,71 -'/Hr 

N. Thus each state 



1/4 ~ 



VI. MQST AND THE NEGATIVE MASS 
QUANTUM PARTICLE 

The equivalence between the repulsive and attractive 
BEC cases, given by Eq. (|20p . means, for instance, that 



9 



the repulsive BEC also contains the Schrodinger cat-state 
(see also Ref. (29f), which is an excited stationary state. 

The double-well potential V — V- (z) for the classical 
phase = TT is responsible for the MQST states of a 
repulsive BEC in a symmetric double-well, predicted in 
Ref. [27] and observed experimentally in Ref. [T^ . There 
are two types of the MQST: the phase-locked (0 « tt) 
states and the running phase states {(j) cx t). Consider 
the symmetric double- well trap (e — 0). Following the 
arguments of Ref. [131 in the classical case, i.e. using 
H = E and Eq. (|23|) one arrives at the equation 



- E 



(43) 



which results in inaccessible z = Q region for \E\ > 1. 
Since for iV 1 the energy satisfies V+(z) < E <V-{z), 
this is possible only when the potential V-{z) = + 
Vl — z'^ has a local minimum at z = (i.e. it is an in- 
verted double- well) and the energy line crosses it. Hence, 
the MQST is possible only for 7 > 1. The mean- field 
condition for the MQST reads 7 > 7c where ^ 



of BEC atoms. This fact makes possible the experimen- 
tal observation of the phase- locked mean- field MQST. A 
similar conclusion was made before using different ap- 
proach [s^, [3l|. We note also that the vr-phase MQST 
was analyzed in Ref. (soj by considering the numerical 
eigenvalue spectrum (see also Ref. [2^) and was related 
to the appearance of the quasi-degenerate energy dou- 
blets. 
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FIG. 5: (Color online) An illustration of the quantum condi- 
tion, i.e. n/hfl- <^ 1, for the phase-locked MQST states for 
several values of the number of BEC atoms A^. 



7c 



,1 + ^1 - Z2(0) COS 



(44) 



One can easily verify that the mean-field critical value 
always satisfies 7c > 1, which is easily seen by rewrit- 
ing Eq. dMl) as 72(0)2/2 = 1 + \/l - z(0)2cos(^(0) and 
noticing that if 7 < 1 the functions on the l.h.s. and 
on the r.h.s. have no intersections for < \z{0)\ < 1. 
Note that for cos 0(0) > we get 7c > 2. For 7 > 2 the 
MQST condition reads |z(0)| > Zc where Zc is the solu- 
tion of Eq. (HH) . On the other hand, for cos (j){0) < and 
1 -I- I sin 0(0) I < 7 < 2 there is an interval of the initial 
population imbalance for MQST: zi < |z(0)| < Z2, where 
Z1.2 solve Eq. 

If repulsive BEC is prepared in one well of the double- 
well trap and the initial phase 0(0) « tt one can use the 
negative mass Schrodinger equation (|18p to explain the 
quantum dynamics. In this case the quantum potential 
in the Fock space V = V-{z) has the double-well form 
with the quasi-degenerate energy levels, which reflect the 
existence of two classical fixed points Zs*' = ±-\/l — 7"^. 
The degeneracy is estimated as twice the tunneling coeffi- 
cient of Eq. ([39| . Therefore, the observation of the mean- 
field MQST with phase (0) « tt is subject to the quantum 
condition that the oscillation time of the effective quan- 
tum particle in one of the wells is much less than the tun- 
neling time between the wells of the double- well potential 
V-{z), i.e. in the energy terms 2k ^ /ifi- — 2h^JaZhZ 
or 



(7^ 



1) 



exp 



N 
'2^ 



^3/2 



1 



N 



(45) 

which, as is demonstrated in Fig. [51 is satisfied for all 
7 just above the critical value, even for small number 



Numerical simulations show that the quantum average 
values defined as [33 



9 ^ 

k=0 



N-1 



(0) = arg(e^P) = arg ^ C^C, 



k=0 



k+1, 

(46) 

initially follow the mean-field dynamics, see Fig. [Hi The 
atomic distribution in the Fock space remains localized in 
one of the wells of the inverse quantum potential — (z). 
Fig. [Zl^a). The distributed quantum phase 0fe, defined 
as 0fc = arg(C^Cfc+i) remains very close to tt, see Fig. 
[It^b). For long times the classical oscillations are subject 
top collapses and revivals, see Ref. The numeri- 

cal method of propagating the Schrodinger equation is 
adopted from Ref. [13] (see also Ref. [Toj). 

To get a quantitative estimate of the valid- 
ity the classical dynamics, one can use the stan- 
dard deviations of the quantum variables, defined as 
Az = 2Ax = 2\\x^{x)\\ = [{{x - (a;))2)]5 and 



Ae'P = \\e'P-{e'P)\\ 



(47) 



(we assume that ](A^, Oj^*)] = 0, in the case of 
1(A^,01^')1 ^ and 1(0, TV]*)] = one can use the op- 
erator e~*P instead with similar result). From Eq. (|47|) 
one concludes that for the average phase defined in Eq. 
(j^S)) we have (e'^) = ae*^'^^ with some a < 1. The quan- 
tum dispersions satisfy the uncertainty relation 



Ae'P 



:Az > h. 



(48) 
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FIG. 6: (Color online) The phase-locked MQST state dy- 
namics: quantum averages (black lines) vs classical dynamics 
(grey lines). The dashed line in the upper panel gives the 
classical stationary point Zs — 0.866. Here = 1000, 7 = 2 
(whereas 7c = 1.23) and e = 0. The initial Gaussian distribu- 
tion of width Az = 0.12 centered at the point zq — 0.78 and 
with phase (j) — n was used. 
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FIG. 7: (Color online) The atomic number distribution (a) 
(represented by vertical bars) of the phase-locked MQST state 
and the distributed quantum phase (b) (given by bars). The 
line in panel (a) gives a schematic portrait of the inverse quan- 
tum potential -^.(1 - 2k /N). 



For the simulations presented in figure [SI the standard 
deviation Az decreases from 0.122 settling to 0.096 while 
Ae'P grows from 0.0479 to 0.268 (the average \ {e^P)\ de- 
creases from 0.999 to 0.963). In this case the l.h.s. in the 
inequality (|48|l grows in the result of evolution by less 
than one order of magnitude as compared to the initial 
value on the order of the r.h.s.. 

It was the running phase MQST which was observed 
in the experiment of Ref. [19]. The classical dynamics 



of the running phase MQST is well understood [27|. We 
have found that the quantum correction to the classical 
running phase is in the form of the quantum collapses 
and revivals, see Figs. [5] and [51 Though we show the 
numerical simulations for small number of atoms N — 
200, the time of occurrence of the first collapse of the 
running phase does not seem to depend on the number 
of atoms but on the initial conditions, for instance the 
phase (the subsequent quantum revivals and collapses do 
depend on N: we have found no revivals for N — 1000 up 
to T = 200). This effect can be observed in an experiment 
if the dynamics is followed for longer times. 
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FIG. 8: The quantum collapses and revivals of the running 
phase MQST state. We use N = 200, 7 = 5, £ = 0. The initial 
condition is a Gaussian with the average 20 = 0.91 (the y-axis 
is cut below the upper limit of z for better visibility of the 
dynamics at large times), 0(0) = (thus 7^ = 3.35) and the 
Gaussian width a = O.IA'^. 



For 0(0) = the first collapse of the running phase 
occurs about the value 0(t) — Stt, which is close to final 
value of the phase in the experiment of Ref. [HI on the 
MQST (this depends on the initial conditions: for 0(0) = 
0.757r the first collapse occurs at 0(t) « IItt, see Fig. [TT|) . 
The growth of the quantum average (0) is interrupted 
by plateaus of constant phase, while the classical phase 
follows the linear growth, see Figs. [TU|and[TT] In the case 
of 0(0) = 0.757r, Figs. [glandfTTl the initial state is more 
classical, i.e. the corresponding uncertainty relation 
is approximately equality at t — (whereas in the case 
of = 0, initially, the l.h.s of Eq. is larger than the 
r.h.s. by an order of magnitude). 

Our principal result is that the first quantum collapse 
is associated with an exponential growth of quantum fiuc- 
tuations of the phase distribution, see Figs. [T^ and fOl 
which reach a maximal value at the first occurrence of the 
quantum collapse. In this case, the l.h.s. of the inequal- 
ity (j48p grows by more than two orders of magnitude 
reaching the value of order one at the first collapse (the 
average |(e*^)| decreases from 0.97 to 0.03). The growing 
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FIG. 9: The same as in Fig i] (iV = 200, 7 = 5, e = 0) 
for the Gaussian initial condition with the average zo ~ 0.8 
and the initial phase (/)(0) = O.TSvr (thus 7^ = 1.8) and width 
a = 0.05Af. 



fluctuations are also seen in the experimental results on 
the running phase MQST presented in Ref. p^ . 
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FIG. 11: The same as in Fig. [TO] but for (/)(0) = O.TStt, i.e. 
corresponding to Fig. [9] 
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FIG. 12: The dispersions of the quantum distributions of z 
and e^^ corresponding to Figs. [8land [T0l The insets show the 
full picture of time dependence of the dispersions. 



FIG. 10: The classical dynamics (thick dashed lines) and the 
corresponding quantum averages (thin solid lines) at the first 
occurrence of quantum collapse of the running phase MQST 
state of Fig. [8] 



Finally, we note that the running phase MQST can be 
described by the WKB approach by setting in Eq. 
= (/)(r) (i.e. the classical solution) and solving the time- 
dependent Schrodinger equation analogous of Ec^. (fT5|) . 
In this approach, it is the running classical phase what 
makes the potential p6p . seen by the effective quantum 
particle, asymmetric in cc, since now we have Vmqst = 
Vq{x, (j)) + <j)x. This gives an elementary explanation why 
BEC is trapped in one well of the double-well trap. 



VII. CONCLUSION 

We have proposed an analytical approach for descrip- 
tion of quantum phenomena in the system of large num- 
ber of interacting identical bosons occupying only few 
modes (two in the present study). The method links 
the many-boson system with the dynamics of a single 
quantum particle in a potential, where the normalized 
occupation numbers of the different modes in the Fock 
space serve as the particle co-ordinates. We have used as 
the example the well-known two-mode model, describing, 
for instance, BEC tunneling in a double-well trap, i.e. 
the boson Josephson effect. The method allowed to ac- 
count for the mean-field stationary points and their sta- 
bility, trace analytically the transformation of the ground 
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FIG. 13: The dispersions of the quantum distributions of z 
and e*'' corresponding to Figs. l9land [TT] 



state of the system of identical bosons in both attractive 
and repulsive case, derive the quantum fluctuations of 
the number of atoms in the ground state, relate the ap- 
pearance of the macroscopic quantum self-trapping phe- 
nomenon to the double- well shape of the potential for the 
effective quantum particle and give a quantum explana- 
tion of the phase-locked and running phase self-trapped 
states of BEC in the double-well trap. We also predict 
a new phenomenon - quantum collapses and revivals of 
classical running phase of the macroscopic quantum self- 
trapped state. 

Our method awaits other important applications, 
where few-mode boson models naturally appear, includ- 
ing the theory of molecular- atomic coherence in BEC 
[34| . the quantum model of nonlinear intraband tunnel- 
ing of BEC in optical lattices 0, [13] and many others. 
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APPENDIX A: THE FULL TWO-MODE 
BOSON-JOSEPHSON MODEL 

One can show that the full two-mode boson model de- 
scribing BEC in a double-well trap can be cast in the 
(dimensional) form 



H 



SVa\ai - [Jiin + Jnoni(A^ - l)](ala2 + 4«i) 



Pi - P2 



y [a[a2 + a^ai 



where some scalar iV-dependent term has been discarded. 
The coefficients are given as Jjin — 5E/2, 



Jnon\ = -g J d^ X(^a(x)(^^(x), 

/ d^^x^Kx), 



Pi ^9 



P2 



(A2) 



(the subscripts a and (3 are permutation of the hst {L, R} 
and the functions give the localized states in the left or 
right well defined by the appropriate linear combinations 
of the ground state and the first exited state, see section 
HI)) . The derivation is similar to that of Refs. [H, [2^ 
and is omitted. However, with the help of numerical 
evaluation, one can verify that for all double-well traps 
with two lower degenerate levels Ei and E2 satisfying 
the inequality 5E = E2 — Ei <^ E^ — E2 the coefficients 
satisfy 

Pi > P2 > iJnonll- (A3) 

The coefficients Jun and Jnoni(A^— 1), however, can be of 
the same order. Discarding the small terms and dividing 
the Hamiltonian by the quantity [Jiin+>^noni(-^~l)]-^ one 
gets the reduced model given by Eq. ([T]) with a different 
definition of the parameters. 



APPENDIX B: PROOF THAT THE LOCAL 
MINIMUM OF -V-{z) AT THE NEGATIVE z IS 
THE ABSOLUTE MINIMUM 

Denote z_ and the negative and positive lo- 
cal minima of —V-{z) for 7 > 7cr. Using 
ez = z^{j — (1 — z^)"^/^) from Eq. (HH) into the expres- 
sion for the quantum potential we get 




-7 I +\/i-4- 



(Bl) 



Using Eq. (pij) again to obtain (1 — z^) > 7 and 

(1 — z^)^^/^ < 7, we arrive at the estimates 



Vl(z+)<^^, V^{z^)> 



1 



where the equality sign is for e 
function /(z) = (1 - zV2)/Vl~I 



(B2) 



1 - zt 



0. Now, since the 
^ is monotonously 



growing for < z < 1, one needs to compare just the 
values |z_| and 2;+. But it is evident from Fig. [TJb) 
that Z4. < \z-\ (the equality for e = 0) which leads to 
V-{z_) > V-{z^). Indeed, z± are the intersections of the 
curve Ihs(z) = 7 — (1 — z^)~^^^ with the two branches 
of e/{2z), while its intersections Zq*'' with the z-axis are 
equal in the absolute value. But, since the derivative of 

Ihs(z) has the sign opposite that of z, we get z_ < Zq -* 

A+) 



(Al) and z_|- < Zq what completes the proof. 
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APPENDIX C: THE DETAILS OF THE 
ANALYSIS OF THE GROUND STATE OF BEG 

To find the ground state of large BEC (TV > 1) one 
needs to compare just the zero-point energies at the ex- 
tremal points (which is the classical energy of the sta- 
tionary point, Eqs. ([24]) and ^E^). Using Eq. we 
get 



74 



1 



(CI) 



is the corresponding extremal point of 



where z± 
V±{z). ^ 

Consider the symmetric case e = and a repulsive 
BEC 7 > 0. The extremal points are z+ — and z_ = 
for 7 < 7cr = 1, while = 1 — 7"^ otherwise. Using 



them into Eq. (ICip we get 



e, the ground state of BEC for e ^ is given by the 
ground state of Eq. pT|) with positive mass. Indeed, we 
estimate from Eq. (jCip using ([2T|l : 




F_(z_)> 



(C4) 



We arrive at the needed inequality V+{z+) < V-{z-) for 
all e, since for a repulsive BEC the absolute minimum of 
— V-{z) is negative z_ < (see Fig. [ija)) while z+ > 
Eq. (in]). On the other hand, by Eqs. (Ell) and (EHl), 
for an attractive BEC z+ is the absolute minimum of the 
double well V+{z) (see FiglTJa)) and is positive, while z- 
is negative. 



ep) 



(^+) = -1, 



v_ 



(rep) 



(Z_) = 



1, 7< 1 
27 



-, 7 > 1 



(C2) 



Thus the ground state of a repulsive BEC in the double- 
well trap is given by the equal distribution of atoms be- 
tween the wells with the zero phase difference. 

On the other hand, for an attractive BEC in a symmet- 
ric double-well trap, using Eq. ()Cip we obtain z- — 
and z+ = for I7I < 1, while z"^ = 1 — 7^^ otherwise. 
Eq. dnH) then gives 



V_ 



(atti-; 



-1, l7l<l 
l7l > 1 



yi""'')(z_) = 1, 



(C3) 

where now 7 < 0. Thus the ground state of an attractive 
BEC also has equal distribution of atoms between the 
wells of the double-well trap (with zero phase difference) 
for I7I < 1. On the other hand, the ground state of an at- 
tractive BEC for I7I > 1, when there are two local bound 
states (in the classical case Zg > and Zg < 0) corre- 
sponding to unequal distributions, is given by the ground 
state ip{z) of the effective particle in the double-well trap 
V+(z), see Eq. (|18p . In this case, the classical stationary 
points feature the spontaneous symmetry breaking. 

Finally, since for arbitrary e > the only change is in 
the position of z± and Eq. (jCip is derived for arbitrary 



APPENDIX D: THE RELATIVE ATOM NUMBER 
FLUCTUATIONS IN THE POSITIVE MASS 
GASE 



We have {{x — Xg)^) = cr+/2, where 



2N 



{l~Zg)'/' 



7 + (1 - z2)-3/2 



1/2 



(Dl) 



From Eq. (|2ip for positive mass case we get (1 — z^)^/^ 
2I 7. Now using Eq. (p9)) we obtain: 



2z, 



-7 = 1 



[1 + 7)^ 



7 + 



2z,, 



7 =(1 + 7) 



8(1 + 7)3 



Therefore 



= — (l+7)-i/2 



l-i!^i±^ + 0(e^) 
16(1 + 7)3 ^ ' 



(D2) 
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